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1 .  INTRODUCTION 


Extensive  research  has  been  done  in  the  past  on  various  problems 
connected  with  the  classical  multivariate  regression  model  since  this 
model  plays  a  very  important  role  in  many  problems  like  prediction.  The 
multivariate  regression  model  is  nothing  but  a  model  with  correlated 
univariate  regression  equations  with  a  common  design  matrix.  But,  there 
are  many  situations  when  it  is  unrealistic  to  assume  that  the  design 
matrices  are  the  same.  One  such  situation  is  when  some  of  the  observations 
on  certain  variables  are  missing.  This  situation  has  been  dealt  with  in  the 
statistical  literature  (e.g.,  see  Srivastava  (1966)  and  Trawinski  (1961))  to 
a  limited  extent.  Another  situation  is  when  the  design  matrices  of 
different  regression  equations  are  not  the  same  but  none  of  the  observa¬ 
tions  are  missing.  For  example,  the  same  independent  variables  may  not 
be  good  to  predict  each  and  every  dependent  variable. 

In  the  sequel,  we  will  refer  to  the  model  based  upon  correlated 
univariate  regression  equations  as  the  correlated  regression  equations 
(CRE)  model.  In  econometric  literature,  the  CRE  model  is  known  as 
seemingly  unrelated  regression  equations  model.  Motivated  by  applications 
in  economics,  Revankar  (1974,  1976),  Srivastava  (1970,  1973),  Zellner 
(1962,  1963)  and  some  other  econometricians  considered  the  problem  of 
estimation  of  parameters  under  the  CRE  model  when  the  underlying 
distribution  is  multivariate  normal.  Recently,  Sarkar  and  Krishnaiah 
(1984)  considered  the  problem  of  estimation  of  parameters  under  the 
CRE  model  when  the  underlying  distribution  is  elliptically  symmetric. 
Approximations  to  the  distributions  of  the  regression  vector  under  the 
CRE  model  were  discussed  in  Maekawa  (1982)  and  Kariya  and  Maekawa 


Kariya,  Fujikoshi  and  Krishnaiah  (1983)  considered  a  model  based 
upon  two  correlated  multivariate  regression  equations  and  they  refer 
to  it  as  the  correlated  multivariate  regression  equations  (CMRE)  model. 
Under  the  above  model,  Kariya,  Fujikoshi  and  Krishnaiah  discussed 
various  procedures  for  testing  for  the  independence  of  the  two  sets 
of  variables  and  also  derived  the  asymptotic  distributions  of  the 
statistics  associated  with  the  above  test  procedures.  But,  no  work 
was  done  so  far  on  tests  for  sphericity  under  the  CMRE  model. 

In  this  paper,  we  discuss  asymptotic  distributions  of  various 
test  statistics  for  sphericity  under  a  CMRE  model.  The  likelihood 
ratio  test  for  sphericity  was  derived  by  Mauchly  (1940)  when  the  under¬ 
lying  distribution  is  a  multivariate  normal  with  unknown  mean  vector. 
Lee,  Krishnaiah  and  Chang  (1977)  approximated  certain  powers  of  the 
likelihood  ratio  test  statistic  for  sphericity  with  Pearson's  type  I 
distribution  and  the  accuracy  of  this  approximation  is  good  for  all 
practical  purposes.  If  we  know  in  advance  about  the  structure  of  the 
covariance  matrix,  we  can  take  advantage  of  this  knowledge  to  propose 
more  efficient  estimates  of  the  location  parameters  and  better  tests 
on  these  location  parameters.  So,  it  is  quite  important  to  investigate 
the  structure  of  the  covariance  matrix  of  the  underlying  distribution 
and  the  independence  of  the  two  regression  equations.  The  results  de¬ 
rived  in  this  paper  are  useful  in  studying  the  robustness  of  the  LRT 
test  for  sphericity  when  the  assumption  of  the  same  design  matrix  is 
violated  under  the  usual  multivariate  regression  model. 
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In  Section  2,  wo  give  some  preliminaries  and  state  the  problems  that 
will  be  investigated  in  this  paper.  Throughout  this  paper,  we  use 
an  estimate  of  the  covariance  matrix  which  is  based  upon  the  residuals 
connected  with  the  regression  equations.  In  Section  3,  an  asymptotic 
expression  is  obtained  for  the  null  distribution  of  the  LRT-like  test 
statistics  for  sphericity.  When  the  design  matrices  of  the  regression 
equations  are  the  same,  the  above  test  statistic  reduces  to  the  LRT  test. 

For  large  samples,  the  asymptotic  distribution  of  the  LRT-like  test  is 
chi-square  and  it  is  the  same  as  the  asymptotic  distribution  of  the  LRT 
test  statistic  for  sphericity  when  the  design  matrices  of  the  regression 
equations  are  the  same.  But,  if  we  take  higher  order  terms,  the  expressions 
for  the  distributions  will  be  different.  In  Section  4,  we  derive  the 
asymptotic  nonnull  distribution  of  the  LRT-like  test  for  sphericity  under 
the  CMRE  model  under  fixed  alternatives.  The  expression  obtained  involves 
normal  density  and  Hermite  polynomials.  The  asymptotic  distribution  of 
the  LRT-like  test  under  local  alternatives  is  given  in  Section  5.  The 
expression  derived  in  this  section  involves  a  linear  combination  of  non¬ 
central  chi-square  variables.  The  results  of  Section  3-5  are  derived 
under  the  assumption  that  the  underlying  distribution  is  multivariate 
normal.  In  Section  6,  we  have  shown  that  the  results  of  earlier  sections 
remain  true  when  the  joint  distribution  of  all  the  observations  is  elliptically 
contoured.  Section  7  is  devoted  to  a  derivation  of  the  moments  of  the 
estimate  of  the  covariance  matrix  when  the  joint  distribution  of  the  ob¬ 
servations  on  each  variable  is  elliptically  contoured  but  we  do  not  assume 
that  the  joint  distribution  of  all  observations  is  elliptically  contoured. 

In  Section  8,  it  is  shown  that  the  asymptotic  null  distribution  of  the 
LRT-like  test  statistic  is  a  linear  combination  of  chi-square  variables 
with  one  degree  of  freedom  when  the  underlying  distribution  is  as  assumed 


in  Section  7. 


2 .  PRELIMINARIES 


Consider  two  correlated  regression  equations 


VX1611+E1 


Y2S*X2622  +  E2 


(2.1) 


where  the  design  matrices  :  n*^,  X^:  n  *  r2  are  known  and  are 
assumed  to  be  of  full  column  rank.  The  matrices  9^:  x  and 

S22:  T2*^2  C^e  Parameters  are  unknown.  We  assume  that  the  rows 
of  E * (E^.E^)  are  distributed  independently  as  multivariate  normal 
with  mean  vector  0  and  covariance  matrix  I,  where 


£  - 


11 

*12 

21 

Z22 

(2.2) 


and  £..  is  of  order  p.  *p,.  An  estimate  of  Z  is  -S  where 
ij  i  j  n 


S-  E'E 


Wl  YiW2] 


•WlYl 


Y2Q2Y2 


(2.3) 


and 


(2.4) 


In  this  paper,  we  are  interested  in  investigating  the  asymptotic 

null  and  nonnull  distributions  of  various  test  statistics  associated 

2  2 

with  testing  the  hypothesis  Hq :  E  ■  o  I  where  o  is  an  unknown  con¬ 
stant. 

We  now  discuss  a  representation  of  S  which  is  used  repeatedly  in 
the  sequel  in  deriving  some  of  the  distributions.  This  representation 
is  due  to  Kariya,  Fujikoshi  and  Krishnaiah  (1983)  and  it  is  given  in 


Lemma  2,1 
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Consider  the  transformation 

(2.5) 

where  is  of  order  (r^-r^)  x  Pj  and  is  of  order  (a-r^)  x  p^.  Also 
Zj,:  (n-r j )  *  n  satisfies  ZjZ..  -  T^/Z^*  -  (^  and  is  chosen  in  the 
following  special  way. 

Let  Zgi  n*  (n-r^)  be  a  matrix  satisfying 


W.  =  Z! Y.  = 

i  ii 


V 


Wo* 


v'y  =  T 

cro  V 


n-r. 


(2.6) 


where  ()q=I-X(X'X)  X',  and  X  =  [Xj»X2]  where  A*  denotes  the  Penrose 
inverse  of  A.  Further,  let  be  the  projection  matrices  onto 
L(X)  nL(Qj),  (j  =  (,2)  where  L(A)  denotes  the  column  space  of  the  ma¬ 
trix  A.  Also,  Jet  Zj  be  a  matrix  satisfying 

Q  =Z  Z:  and  Z'Z  =  I  .  (2.7) 

J  J  J  J  J 


Tlien  choose 

Zl  =  (Zr  ZQ),  Z2  -  (Z2,  ZQ).  (2.8) 

It  is  easy  to  verify  that 

Z'Z  =1  for  n.  =  n-r.  and  Z  Z.'  -  Q.  , 
i  i  n.  '  1  ill 

where  is  given  in  (2.4).  Note  that  under  H^,  the  rows  of  are 

independently  and  identically  distributed  (i.i.d.)  as  multivariate  nor- 

2 

mal  with  mean  vector  0  and  dispersion  matrix  a  I  .  Also,  W.  and  W. 

Pi  1  2 

are  distributed  independent  of  each  other. 

From  (2.3)  and  (2.4)  we  obtain 
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I 

I 

a 


a 


r 

I- 


i- 


where 


U' w 
wlwl 


wiz*z2w2 


W2Z2Z1W1  W2W2 


G  +  B, 


(2.9) 


ru« 


<U1  U2> 


(2.10) 


MiM! 


MiKM2 


M2M2 


,  K“ziz2- 


(2.11) 


These  results  can  be  summarized  in  the  following  lemma: 


Lemma  2.1.  The  matrix  S  defined  by  (2.3)  can  be  written  as  G+B,  where 

G  and  B  are  defined  by  (2.10)  and  (2.11)  respectively.  Under  Hq,  G  is 

2 

distributed  as  wp(ng>  °  Ip)>  where  p  *  P1+P2  and  ng  *  n-rQ  and  row^M-[^ 

2 

is  distributed  as  N.  .  [0,1  ft  a  I  ],i*l,2,...  where  M 

(r0  i,Pi  ~  r0"ri  Pi  1 

is  defined  by  (2.5).  Also,  G  and  B  are  distributed  independent  of 

each  other  under  Hq. 

3.  ASYMPTOTIC  NULL  DISTRIBUTION  OF  THE  LRT-LIKE  TEST  STATISTIC 
The  hypothesis  Hq  can  be  tested  by  using  the  statistic 


ls|n/2 

(trS/p)"^2  ' 


(3.1) 


When  X^”X2,  the  above  statistic  is  the  likelihood  ratio  test  statistic 
for  sphericity.  We  will  derive  the  null  distribution  of  T«a.T  where 


T  * -2  log  A 

»  n[p  log  tr  S  -  log  |  S  |  -  p  log  p] 


(3.2) 
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and  a=(nQ/n)£l,  under  the  assumption  that  K=Z|Z230(1)  as  n^-*-' 

Here  we  note  that  A  is  not  the  LRT  test  statistic  when  X  +X  . 

1  '  2 


so  that 


V  =  *^0  ^n~  ~  v*1  * 
0  n0  P 


(3.3) 


0  •  V°\ 

/no 


(3.4) 


S  -  n  a  (I+A), 


(3.5) 


where 


A  -  (- 


log  | S |  =  p  log  (nQa  )  +  log  |l+A| 

logtrS  -  log  (nQa2)  +logp+log  (l +-ElA) 

P 

From  (3.2),  (3.6)  and  (3.7)  we  obtain 


T  “  TQ  +  —~2  (Tx+T  ), 
/n0 

where 


I0-^Ur 

r 


3 

tr  VJ] 


T  =  — -  f  (Cr  ^ 
1  1  2 
3o  p 


V"7  (tr(VB)  -  ~--tr  B  ] , 
o  F 


(3.6) 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


(3.11) 


The  characteristic  function  of  T  Is 


E[e1C  TJ 

E[e  0  {1  t  }]  +E[e  °  ~  +  CKn^1) 

*^0 


( t )  +  <j>2(t)  +  0(nQ  ) 


(3.12) 


where 


it  T 

^(t)  =  E(e  °  {1+—  T^] 
it  T 

^2(t)  *  E[e  °  —  T2]  . 


(3.13) 


(3.14) 


Note  that  the  characteristic  function  of - -  log  [  — lfi-L  ]n^ 

-1  2  "  (-^)P 
is  (|>  (t)+0(n  )  where  G~W  (n  ,o  I  ).  So  we  know  v  p  1 

1  U  pup  r 


^(t)  =  (l-2it)"f/2+0(n"1) 


(3.15) 


where  f  *  (p(p+l)/2)  -  1. 

Next  we  consider  <{> 2  ( t ) .  Taking  expectations  with  respect  to 
only  yields 


E[trVB]  =  E(tr(VuM^M1+V12M^K'M1+V21M^KM2+V22M^M2)] 


tr[V11E(M[M1)  +V22  E(M^M2)] 
a2(r0-ri)tr  Vu+o2(r0-r2)  t  rV^ 


(3.16) 


E[tr  B]  =  E[tr  MjM^  +  tr  M^] 

*tr[a2(r0-r1)Ip^+a2(r0-r2)Ip2] 
-  [o2(r0-r1)p1+a2(r0-r  )p2]. 


(3.17) 
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From  (3.14) 


it  T. 


<1>2  <t)  "  £le 


where 


it 

(r2_ri 

P  Q2 

h 

V,  „ 

11 

12 

V 

1  21 

22 

:  V  pi”pj 


Now,  note  that  the  limiting  distribution  of  V“(v^)  in  (3.3) 

_  _  4_4 

distribution  of  V»(v^j)»  where  vi^~N(0,2o  ),  v^~N(0,o  )  i 

vi.  (i  < j)'s  are  all  independent,  and  that  the  density  f  of  V 
expressed  as 


f(v)  =  f0(v)  +—f1(v)+^-  f2(V)  +  ... 
\  0 

where  fQ(V)  is  the  p.d.f.  of  V.  Next,  let 


ixp 


v*'  =  (v*,v*, . . .v*2) , 


(VU p v22 1  ‘  *  VPP  *  V12  p '  *  *  Vlp  ’ V21 '  •  *  *  V2p '  * '  *  Vpl 


* 

Then  the  limiting  distribution  of  v  is 

N  2  (0,  Da) 

P 

A  A  A  A 

where  Do*Diag(2a  ,...,2o  ;  a  ,...,a  ). 

Further  note  that 


*v  P 

pp-1 


vij+F{(p-1)^vu-.Lv 


2o  i*j 


i*j 


ii  3j 


)] 


1  *’  *  * 

=  — t  v  A  v  ,  say 

2o  ' 


(3.18) 


is  the 
jf  j  and 
can  be 


(3.19) 


(3.20) 


(3.21) 


where  A*  depends  on  p. 
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So, 


p1 

It  !  1  **  -1  ft  *  ,  r*  F  *  *  -1 

♦2(t )  -  c  oxpl-  i  V  (d/--1-;-  A)v*Hp2  [  vi-pi  I  vi]dv  +0(n  l) 

✓n7  o'  p^+1  D.+l 


Pl+1 


-  0  +  o(n^) 


(3.22) 


Using  (3.12),  (3.15)  and  (3.22),  we  get 


>(t)  =  (l-2it)“f/2+0(n"1) 


(3.23) 


where  f  «  (p(p+l)/2)  -  1.  Mow,  inverting  the  right  side  of  (3.23) 
yields  the  following  expression  for  the  asymptotic  distribution  of  T: 

Pr  (T<x)  -  Pr  (Xfix)  +  (KnjJ1), 
where  f  *  (p(p+L)/2)  -  1 

4.  ASYMPTOTIC  NON-NULL  DISTRIBUTION  OF  THE  LRT-LIKE  TEST  STATISTIC 
UNDER  FIXED  ALTERNATIVE 

Let  us  consider  the  alternative  not  H^.  Since  the  test  sta¬ 
tistic  A  is  a  function  of  the  eigenvalues  of  S,  we  can  assume,  with¬ 
out  loss  of  generality,  that 


E  “  D^  =  Diag(l^ , ^2 , • • • , Ap)  ■ 


(4.1) 


Also,  G,  B,  M^'s,  Oq  are  defined  as  in  Section  2.2. 
Let 


v  ■  ^  V 


(4.2) 


Under  H,  ,  G~W  (n„,D.).  Now 

1  p  U  A 


S  -  n  (D  +  -'L  +  J-  ) 
0  X  ^  n0 

"  W1  +  A)’ 


(4.3) 


where 


A  =  <—  D~XV  +  ~  D~LB) 
/nQ  0 


log  |  S I  »  p  log  nQ  +  ^  log  +  log  1 1  +  A| 


i=l 


tr  S  =  nQ(  i  A.+trC),  C  = 

i=l  \  0 


=  n.pX(l  +  , 

PX 


where  pX  =  X,-r. .  .+X  .  Also, 
1  P 


p  log  tr  S  =  p  log  rig  +  p  log  p  +  p  log  X  +  p  log  (1  + 


Hence 


T  *  n[log 


^ - +P  log  (1+-^^)  -log|t+A|]. 

h  x  px 

l-x  1 


Now,  let 


^  [iT-l0gF7] 
A1 


=  </n^  [p  log  (1  )  ~  log  |l+A|] 

pX 


T0  *j=.  <VV  +  ^“O1’' 
% 


where 
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?«i- 


Ta  •£«=  - 


The  characteristic  function  of  T  is 


(4.10) 


where 


4>(t)  =  4/  (t)  +  ip»(t)  +  0(n  1) 

1  4  O 


<l<1(t)  -  E[e  0  (1  +  —  Tj)] 


(4.11) 


itT0  it 

*2(t)  =  E[e  T2]. 

0 

Defining  v*  as  in  (3.19),  we  see  that  the  density  of  v*  i8 

p2xl 

N  2(0,A)  +  0  (nQ^2),  where 
P 

2  2 

A  ■  Diag(2A  , . . . ,2A  ,A  1  , . . .X, A  ,A  A  .....A  A  ). 

i  pl2  lp21  p  p-1 


a,  -  (=  "  T~)  .  i  *  l(l)p 
1  A  i 


(4.12) 


(4.13) 


(4.14) 


(4.15) 


n  ~  ( a  ^ ,  a  2 , . . . ,  a^  ,  0 , 0 , . . . ,  0 )  • 


(4.16) 


From  (4.8),  we  know  that  TQ  =  a'v*.  Also  from  (2.9),  ^  can  be  written 

2  2 

as  where  elements  of  fl:  p  *p  depend  on  A,,..., A  and  p.  Then 

1  P 

from  (4.12), 


ita'v* 

4»  (t)  -  E[e  (1  +  —  v*'Qv*)] 

t2  ,  S' 

- 7-  a'Aa 

L  •*  — 


3  3 

11+—  tr  <0A)  +  -^~  a’AnAa+-^ - 3 - - -  ]  (4.17) 


r/>- '  ■  '  ■  ’ 

S 

L* 

13 

[ 

We  now  consider  ip2(t).  Taking  conditional  expectation  with  re- 

B 

spect  to  M^'s  only  yields 

►  * 

P1  P 

w 

VV  ‘  E|  +,  i  *ibil> 

1=1  i=p^+l 

*  . 

1 

r1  E 

“ — ; 

=  ^  ai(Vri)Xi  +  i  ai(Vr2)Xi 

i=l  i=Px+l 

Jl  E 

1 

"  i  aiXi  +  (ro“r2)  4  aixi 

i=l  i=p^+l 

*• 

=  K1(r0,r1,r2,ai' s,Xi's,p1,p2) 

i 

5  1^,  say  . 

(4.18) 

Hence 

i 

,  v  it  „  _r  it  a'  v*, 

^2(t)  =  Ki  E^e  ~  '  1 

^  t2 
,lt_K  e-2a.’Aa 

1 

(4.19) 

- 

i 

0 

Finally,  from  (4,11),  (4.17)  and  (4.19).  we  have 

-  1 

>• 

2  2 

-  -LJ-  .  3 

*(t)-e  2  I1+-:  g1+^-  g3]  +0(n"1) 

(4.20) 

where 
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Note  that  under  H^,  i  t  0,  so  that  inverting  the  rightside  of 
(4.20),  we  have  the  following  theorem  for  the  asymptotic  distribution 


of  T. 

_P 

;  , _  t  ^ 

Theorem  4.1.  The  distribution  function  of  T*  ■  T/t  -  /n-  (—  T  -  log trr— ) /t 
"  u  n  ii 

under  can  be  expanded  for  l.arge  n  as 

(3) 

Pr(T*  •  x]  =  >l>(x) - -  [g1  (x)/t  +  g3  —  3^  ?  +  °(nQ1) 

where  $^(x)  is  the  jCh  derivative  of  the  standard  normal  distribution 
function  <b (x)  ;  and  g  ,g  and  t  are  given  by  (4.21). 


5.  ASYMPTOTIC  NON-NULL  DISTRIBUTION  OF  THE  LRT-LIKE  TEST  UNDER 
LOCAL  ALTERNATIVES 

We  assume  the  same  structure  of  Y.  as  in  (4.1),  but  we  consider 
local  alternatives 


v  i  ■  i(i)p 

*\> 

where  (J.'s  are  not  all  equal.  Thus,  under  H„ , 

l  0 


(5.1) 


D  -  XI  +  — 

A  p  /— 

% 

where  D  «  Diag(0  ,0  ,...,0  ).  Under  H  , 

i  i  p  u 


C~W  >n,<Xl  +  ~ )]. 

p  U  p  ) — 

s 

Define  V  as  in  (4.2),  where  D^  is  given  by  (5.2).  We  then  have 


(5.3) 


S  =  nnX(I  +A) , 

0  p 


(5.4) 


where  A  ■  ((V+D(J)/A  i/n^)  +  (B/A  Oq)  .  Expanding  log  |s|  and  log  tr  S 


in  the  same  way  as  before,  we  get 


T  *  -  2a  log  A,  a* —  <  1 
n 


where 


T0  +  — (T1+T2) , 
/n0 


(tr  V  +  I  e  r 

T0  =  ^2  ftr(V+D0)2 - i - ] 


'  +  5  e.: 


T  =»— =—  [- 
1  3  1 

3AJ 


- tr(V+DQ)3] 


tr  B(tr 


T2  — ~  [tr  B(V+D0)  - 


V+  fo  ) 
1  * 


The  characteristic  function  of  T,  as  before,  can  be  written  as 


<t>(t)  =  *1(t)  +*2(t)  +0(nQ1), 


where  ^(t),  4>2 ( t)  have  the  same  expressions  as  in  (3.13),  (3.14) 
respectively. 

Now,  ^(t)  +  0(ng3)  is  the  characteristic  function  of 

-2  n  /n  log [  --J^r  -']n^2,  where  under  H  ,  G~W  [n.,(Al  +  -33— )],  and 
u  tCT  \ P  B  p  0  p  r— 

(  P  )  /no 

hence  <^(t)  is  known  as  (see  Fujikoshi,  1981) 

2 

<*>!  (t)  -  Mt;A2/A2)[H-  ~  l  A"3b.(l-2it)"j)+0(n"1)  (5. 

j=0  J 

2 

where  ip ^ ( t ; A)  is  the  characteristic  function  of  a  noncentral  x  var 

2  2 

iable  with  f  d.f.  and  noncentrality  parameter  6  /A  ,  and 


21  2-1  2 

°  D0  -  P  (tr  V  3 


bn  “  Xl2  tr  ~  3P  1(tr  Dn)tr  D2  +p'2(tr  DQ)3] 


and  A  is  given  by  (4.14)  and  (5.1). 


Since  A^  =  A  +  0( -  ), 

A  DQ  =  Diag(2X2,...,2X2;X2 . X2) 


(5.16) 


The  limiting  distribution  of  V  is  that  of  Z=(z^j)  where  Z^~N(G^,2A  ) 
2 

Z^~N(0,A  ),  i  t  j  and  Z^(i<J)'s  are  all  independently  distributed. 
Hence,  for  large  n 


v*~N  2(V  Dq), 


(5.17) 


where  and  are  given  by  (5.15)  and  (5.16)  respectively. 


Hence 


<j>~(t)  - 


.  (r„-r  )  ,t 

- - I-  c'v*  dv* 


(^T)P2|d0|1/2 


where 


2  *  <P2  ~p  ,_P1 -p0  ’  0  2  > 

l*p^  1 1  v2  p  -p 


(5.18) 


( r2_ rl )  it  i,  it  „  .  i-l/2  r  1,  ,_-l  .„-l  ... 

-  —  U--2  Vo1  expt'  2(1J~0DO  -0_~G^  V1?  Yg 

/no  A 


(5.19) 


where 


o  /n“l  it  .  . -1 

n  -  (Do  -  TT  V 

A 

v  =  hd"1^ 

-G  0  -.6 


(5.20) 


l  (D„Ajr  u. 


(5.21) 
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6.  THE  DISTRIBUTION  OK  THE  LRT-LIKE  TEST  STATISTIC  WHEN  THE  JOINT 

DISTRIBUTION  OF  OBSERVATIONS  IS  ELLIPTICALLY  CONTOURED 

We  will  first  discuss  briefly  elllptically  contoured  distributions. 

If  the  random  vector  x :  p  *  1  has  the  characteristic  function  of  the 
form  exp(i  t'  p)  ij)  (t1  t)  where  p  and  t  are  of  order  p  x  1,  then  x  is  said 
to  be  distributed  as  eiliptically  contoured  distribution  and  is  denoted  by 
EC(u.E.4>).  Various  properties  of  eiliptically  contoured  distributions  are 
discussed  in  Anderson  and  Kang  (1982),  Carcbanis,  Huang  and  Simons  (1981) 
and  Kelker  (1970).  Now,  let 


where  «*'  -  (vecE')'  -  (e'(1)  .. . .  «e  assume  that 

e*  "  Enp(0,  I^i  I*;  4>)  and  T*  is  proportional  to  E  given  by  (4.1). 

Then,  it  is  known  (see  Cambanis ,  Huang  and  Sitnons(1981)  and 


Anderson  and  Fang  (1982))  that 


d 

E  =  R  U  A, 

where  A' A  ■  £*,  A  :  P  *  P,  U  :  n  x  p,  Vec.  U  -  u^np\  distribution 
function  of  R  Is  related  to  <j<  and  R  Is  independent  of  U  ;  here  u^np^  Is 
np-dlmenslonal  column  vector  which  has  uniform  distribution  on 
the  unit  sphere.  In  addition,  "X  ■  Y"  denotes  that  the 
distribution  of  X  is  the  same. as  that  of  Y.  Let  us 


d  d 

write  A  *  (Aj  A?) ,  P  “  Pj  +  |>2.  Then  we  have  *  R  U  A^ ,  ■  R  U  Aj. 


nxPj  nxp2 


Since 


E’l‘hEl  EiW2N 
K2<V?1K1  E2<>2E2 


(6.1) 


we  get 


and 


trS  -  trCK'jQjEj  +  E^Q2E2) 

-  R2tr(A*jU'QjUAj  +  A£U'Q2UA2) 


I S |  -  lEjQjE,  |  |E!,Q2E2  -  K^Q2Q , E ,  ( E j Q ,  E j )“ 1  E'jQ , Q2E2 1 

-  r2p|a'1u,q1uai  |  |aJu'q2ua2  -  a^u'q^jUAj (AjU'QjUAj)-1  a|u,q1q2ua 


(6.2) 


2' 

(6.3) 


Hence 

igi  d  |A*U*QlUA1|  |A^UaQ2UA2  -  A^U'C^Q^  (AjU'QjUA^'1  AjU'QjQ^  | 

(trS)P  [tr(A;u'i).UAt  +  A’U’Q2UA?)]P 

111  i  (6.4) 

2 

Substituting  this  we  see-  that  A  is  independent  of  R  .  Hence  the  distri¬ 
bution  of  A  will  he  the  same  as  in  the  normal  case.  Thus  the  asymptotic 
null  and  nonnull  distribution  of  A  under  assumption  (6.1)  is  the  same  as 


in  earlier  suctions 


20 

7 .  MOMENTS  OF  THE  ESTIMATE  OF  THE  COVARIANCE  MATRIX  WHEN  OBSERVATIONS  ON 
EACH  VARIABLE  ARE  ELLIPTICALLY  CONTOURED 

In  this  section,  we  assume,  instead  of  (6.1),  that  the  covariance  matrix 

of  e'  is  E  ft  I  where  E  is  given  by  (4.1),  e'  -  (e),...,e')  and 
—  -  •»*  «*P 


e,  -  ECn(0,  \*In;  *)  (7.1) 

where  X*  «  X  ^  ,  j  -  l,...p. 

It  can  be  verified  easily  that 

E(e^Ae1)  -  -24>*  (0)  X*trA  (7.2) 

E(e'Ae  )  -  0,  i  +  j  (7.3) 

Var(e^Aei)  -  8  *' (0)2X*2tr  A2  +  12(t"(0)  -  *'(0)2)X*2  f  a^  (7.4) 

for  A  •  A'  “  (a^) 

Varfe^Aej)  -  4  X*X*U"(0)I  +  ■>'  (0) 2£^a±j 2 1  •  1  +  J  •  <7-5> 

Note  that  S  *  E'E 


Note  that  S  *  E'E 


where 


nt  -  n  -  rt,  n2  -  n  -  r2 
E(s  )  *  0,  i  +  j. 


(7.7) 


! 
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84>’(0)2A*2n  +  12(*"(b)*)A*4  l  qJJ,Z. 

1  i=l 


Var( s  . ) 
11 


"(6)2)a*2  v  ~(1)2 

i 

i  -  1  >  2  #  •  *p. 


J  *-»-»••• 

)  8<|)'(0)2X*2n2  +  12 (<f>*’ vO)  — <|> '  (0)  )2X*2  J  qj2)2, 

V. _  i  =  Pj+1, . . .  ,P 


where  Q.  = 


<»>.  q2  -  «,«>> 


4X*A*[r(0)Eq^)2  +  ♦,(0)4tH  <« 

J  a 

i,j  -1,2 . Px;  1  J* j 


2rr  (1)2 


Varts^j) 


4A*A*U"(0)Iq£;)2  +  +'(0)2n  q^)2i 

1  i  a  UOt  a?B  ali 

i,j  *  Pj*1 . p;  i  /  J 

«x*xn*"(o)^«)2  +  ,‘f 2] 

J  a  a^B 


i  -  1 ,2, . . . ,p^ ;  J  —  p » •  • » » p 
or  i-p  +l,...,p;  j-l,2,...,p^ 


(31 

where  =  (q^  ) 


C  4X*X*(<Ji"(0)  -  4>'(0)2)[q 
\  1  J  a  ' 


(1)2 

aa 


i,j  -  1,2,. . .,P1;  i  /  j 
,2.r  (2)2 


CoV<W 


4X*X*(<fr"(0)-<^,(0)Z)^a 

■J  n 


i , j  -  Pj+l,.  • .  ,p;  i  i  j 


4X*X*(4»"(0)-<fr*(u)2,[q^)q^) 

-*  n 


1  -  1,2,... .p^,  j  * p1+l, . . . ,p. 


(7.8) 


(7-9) 


(7.10) 


ft 
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For  simplicity  of  notation,  let 


r  => 


nl  =  n2  =  n0» 


say 


where  nQ  =  n-r. 
Hence 


E(  )  s  Diag(A  ,  ■  •  » »A  ) 
n0  1  P 

where  =  — 2<J> *  (0) X* ,  i  =  l,2,...,p. 

Let  us  define 


Z  * 


V- 


Then 


Var(Z 


ii 


Var(zij) 


E(Z)  =  0 
pxp 

8^’(0)2A*2  +  12(<T(0)-*’(0)2)A*2[  q^)2/nQ 

i  3  •  •  • » P-j^ 

8<j>'  (0)2X*2  +  12(0"(O)-^(O)2)X*2  y  q<f2/n0 
i  =  Pj+1. . • • ,P 

4A*X*[*"(0)I  q^>2/n 0  +  ^(0)2n  qaf2^1 
1  j  q  aa  u  a?«B 

i.j  -  i(i)  Pl ;  ±t  j 

«;x*l*''(0)Sq<f2/„n'(0)2n  ,<“2/n0] 

J  U  01*  P 

i.j  “  Pj+1,  .  •  •  .p;  i  t  j 

4X*X*[r(0)I  q^)2/n0  +  ♦'(O)^  q^/n,,] 

J  a  a^B 

i  =*  KDPj!  j  =  Pj+1, .  •  •  ,P 
or  is  p^+1 . .  j  *  KDp^ 


(  7.11) 


(7.12) 


(7.13) 


(7.14) 


(7.15) 


r  4 A*AJ (°) “♦ *  c°)  ) a  C  /no 


C°V(Zii,Zjj)  = 


i,j  = 1,2, . . . ,p  ;  1 +  j 


4A*A*(^”<O)-0*  <0)^Z  ^)2/n0 


i,  j  =  Pj+1, .  • .  ,p;  i  /  j 


4»*»}(*-<0)-*-(0  )hl 


i*  l(l)p  ;  j  -Pj^+l,...,?. 


All  ocher  elemencs  of  Z  are  uncorrelated.  Now  make  the  following 

assumptions  on  the  design  matrices  Q  .Q^*. 

Each  of  y  q(1)q<2)/n  ,  J  q<j)2/n  ,  J- 1,2,3  and  H  q(J)2 
‘  ,  ‘cm  au  0  •.  au  0  aB 

a=l  a=l  ar& 

j*  1,2,3  are  of  0(1)  and  we  write  for  large  n 


J  "aa’V^’-  J'1-2’3 


'rtS  '  k2J>  •  3  " 1,2,3 

a  ?B 


y  q(1)q(2)/n  =  K  . 

L  ^aa  ^aa  0  3 


Note  that  the  limiting  distribution  of  Z  is  the  same  as  that  of 
Z  =  (z^) »  where 


z..  ~N(0,2A>), 


N(0, i  +  j 


( 
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r 


and 


1  +  |(<h*,(0)/4.,(0)2-l)K^1) 

i —  1*2,.*. 

1  +  f(*"(0)/<J.'(0)2-l)K^2) 

i  =*  Pj+1,. . .  ,p 


(7.20) 


4>  * 


Also, 


.  »"(0) 

2 


^>1. 


i,j  =  l(l)p1 


,  4"(Q)_ 
<J>'  (o)2 


i.j 


p^+1 f  *  * • tP 


(7.21) 


.  ♦"(0) 
♦'(O)2 


K<3>], 


i » l(l)p,  j  -  px+l, 
J  -l(l)Pl,  i-Pl+l 


•  .P 


•  •  »P  • 


C0,('lf'jj) 

•Vj 

-G-  1-©^ 

/-S 

o  o 

W  V-/ 

hO 

-  1)C,  i/j 

(7 

.22) 

r 

K(l) 

K1 

\ 

i 

,  i.j 

-  KDPl 

c  =  < 

) 

„(2) 

K1 

.  i.j 

■  Pj+1 . . 

(7. 

.23) 

( 

V 

K3 

,  i  * 

l (l) p1 ;  j  =  p1+i,... ,p. 

All  other  elements  of  Z  are  uncorrelated. 

It  is  known  that,  all  fourth  order  cumulants,  if  they  exist,  can 
be  expressed  as  a  function  of  a  single  parameter,  k,  which  characterizes 


the  kurtosis  of  the  distribution. 


For  E  as  in  (4.1),  the  only  nonzero  fourth  order  cumulants  of  the 


elements  of  E  are 


K4  “  E<eti)  “  3  Xl  *  -  HDn 

-  3-\^(^-W'2  -  1),  1  -  1(1)P 

1  ♦'  (O)-4 

k2J2  ‘  E<«u*l|>  -  E(e«1)E<«2J) 

ij  =*  i(i)p,  i  +  j. 

1  2 

All  other  fourth  order  cumulants  of  e^'s  vanish.  Note  that  K^/X^, 

Is  the  kurtosis  of  the  marginal  distribution  of  iC^  component,  and 

i  2 

for  convenience  define  K^/A^  -  3k,  i  *>  1  (L)p;  where  k  characterizes 

the  kurtosis  of  the  distribution.  It  is  clear  from  above  that 

k  ■  -  1),  so  that 

*'(or 


The  parameters  in  the  asymptotic  distribution  of  Z  can  be  expressed 
as  functions  of  k  and  A^'s  only. 

8.  ASYMPTOTIC  NULL  DISTRIBUTION  OF  THE  LRT-LIKE  TEST  STATISTIC  WHEN  THE 
OBSERVATIONS  ON  EACH  VARIABLE  ARE  ELLIPTICALLY  CONTOURED 


Let  SQ  *  s/ng>  then  from  (7.12),  we  have 


S 


0 


D.  +  Z//n 


(8.1) 


1  O'  n/2 

The  LRT-like  test  statistic  for  sphericity  is  A  ■  [ - ]  '  . 

(tr  SQ/p)P  . 

We  are  interested  in  the  asymptotic  distribution  of  T  *  -2  log  A,  Sup- 

2 

pose,  under  Hq,  Aj»...«A  «u  .  Then 


Hence , 


SQ  *  a  ( I+A) ,  A  = 


log | SQ |  -  p log  a  +  log | I+A | 

2  +  »•  A 

log  tr  S.  ®  log  a  +  log  p  +  log  (1  +  — )  , 

u  p 


T  =  n[p  log(l  +  ££•£-)  -  log  |  I+A |  ] , 

P 

n0  n0 

let  i'  =  —  T,  —  <  1,  being  the  correction  factor.  Then  we  have 
n  n  — 


T  =  -L  [tr  Z2  -  Z^~ ]  +  0(n~1/2). 


So  the  characteristic  function  of  T  is 


(8.2) 


<f>(t)  =  Efolt:  T] 


^  -1/2 
E[e  °]+0(nn1/2). 


(8.3) 


where 


To  *  “T  z2 

2o 


(tr  ZY 


(8.4) 


Let  us  write  this  as  z  'A^z ,  where 

~  ,=  <Zll,'”,Zpp,Z12,'">Zlp;Z21 . z2p’  *  *  *  *  zpl . zp(p-l)>  (8,5) 

lxp4 

2 

and  Ag  is  a  matrix  whose  elements  depend  on  p  and  a  .  We  assume  that 
the  errors  are  distributed  as  in  Section  7.  Then,  asymptotically. 


z~N  2(o,  a), 
P 

where  H  =  (w^) 


(8.6) 


2 

wij,s  are  given  in  Section  7,  with  the  restriction  that  *  a  for 
i  ■  l(l)p.  Hence, 


Inverting  the  rightsidc  of  (H.7),  we  get 

Theo rew  8 . 1  _  The  limiting  null  distribution  as  T  is  that  of 

a  linear  combination  of  chi-squares  with  one  degree  of  freedom  and 

the  coefficients  depending  on  the  fourth  order  moments  of  the  obser- 

2 

vations  of  the  parent  population,  which  are  functions  of  o  and  k. 
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